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ABSTRACT 

We analyse the model of stochastic re-acceleration of electrons, which are 
emitted by supernova remnants (SNRs) in the Galactic Disk and propagate then 
into the Galactic halo, in order to explain the origin on nonthermal (radio and 
gamma-ray) emission from the Fermi Bubbles (FB). We assume that the energy 
for re-acceleration in the halo is supplied by shocks generated by processes of star 
accretion onto the central black hole. Numerical simulations show that regions 
with strong turbulence (places for electron re-acceleration) are located high up in 
the Galactic Halo about several kpc above the disk. The energy of SNR electrons 
that reach these regions does not exceed several GeV because of synchrotron and 
inverse Gompton energy losses. At appropriate parameters of re-acceleration 
these electrons can be re-accelerated up to the energy 10^^ eV which explains in 
this model the origin of the observed radio and gamma-ray emission from the 
FB. However although the model gamma-ray spectrum is consistent with the 
Fermi results, the model radio spectrum is steeper than the observed by WMAP 
and Planck. If adiabatic losses due to plasma outflow from the Galactic central 
regions are taken into account, then the re-acceleration model nicely reproduces 
the Planck datapoints. 


1. Introduction 


Recently Fermi has discovered two giant gamma-ray Bubbles (F Bs) that exte nd 


nearly 10 kpc in diameter north and south of the Galactic center (cf . 


Sn et al 


201(1. for more recent analyses, see 


Ackermann et al. 


Hooper fc Slatyerl (2013); 


Dobler et a 


Yang et ah. 


2010 


(1201 4h 


(120141 )). These gamma-ray Bubbles also correlate with the earlier 
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discovered so-cal 


Finkbeinerl f 2004j) and 


ed ’’ microwave haze” observed by the WMAP telescope as described by 


Dobler fc Finkbeinerl fl2008f) . and with the large scale X-ray emission 


region hrst evidenced by analysing the ROSAT 1.5 keV data , which clearly showed the 


chara cteristic of a bipolar flow (see e.g., 


Snowden et ah 


1997 : 


Bland-Hawthorn fc Cohen 


20031 ). A number of models was suggested to explain the 


due to protons-proto n collisions (hadronic model, see e.g., 


origin of the FBs el 

ler 

Crocker & Aharonian 

2011 


Crocker et al. 


2014alJbl) o r due to inver se Compton scattering of relativistic electrons 


(leptonic model, see e.g.. 


Su et al. 


20101 ). 


gamma-ray emission has a cut-off at ~ 100 GeV (see 


from observat 

ions. First, 

the 

Su et al. 

2010; 

Yang et al.. 

2014) 


Then, in the case of inverse Compton origin of the FB gamma-rays electrons have to be 
accelerated there up to the energy 


E^ax ^ ~ 5 X 10^^ eV, 


( 1 ) 


where e ~ 10 ^eV is the energy of the microwave photons. 


Secondly, because of a very shor t lifeti me of electrons, they have to be in-situ generated 


in regions of emission. 


Cheng et al. 


FB edge by shocks while in the model of 


(1201 ill assumed that electron s are accelerated at the 


Mertsch fc Sarkarl (1201 ll) it was assumed that this 


emission is generated by electrons accelerated in the Galactic halo by an MHD-turbulence 
which is excited by a shock propagating into the halo. The energy of this shock cascades 
into turbulence by different processes of plasma instabilities. Interaction of electrons 
with this turbulence leads to stochastic or second-order Fermi acceleration. Alternatively 
particles can be stochastically accelerated in the FBs by interaction with a s upersonic 


turbu lence (shocks) which is excited by tidal processes in the Galactic Center (jCheng et al. 
291211 . 


These models were investigated in the test particle approximation when feedback 
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reaction of accelerated particles on the acceleration mechanism is ignored. Therefore, the 
nnmber of accelerated particles is usually a free parameter of the models. However, in some 
cases when the sources of accelerated electrons are known this number can be estimated 
from kinetic equations that gives additional model restrictions. In the case of FBs there are 
three evident sources of electrons. The electrons can be supplied by: a) Coulomb collisions 
from the FB background plasma, b) by pp collisions in the halo (secondary e l ectrons'), and c ) 


SN Rs in the Galactic Disk. We d iscussed models a) and b) in 


and 


Cheng et al 


(12015. paper ij ). respectively. 


Cheng et al 


( 2014. paper j) 


In 


Cheng et al. 


(12014. paper ll) . we analysed a model of stochastic acceleration of 


electrons from the background p 
effect of plasma overheating (see 


asma and showed that the problem of the model is the 


Chernyshov et al. 


2OI2I ). However, for a specihed set of 


acceleration parameters the in-situ stochastic acceleration is able in principle to provide 
high energy electrons needed for the observed radio and gamma-ray emissions from the FBs. 


In 


Cheng et al. 


(12015. paper Ill) we analysed the hadronic model of gamma-ray 


emission from the FBs when gamma-rays are produced by pp collisions while the radio flux 
is generated by secondary electrons. Owing to low gas density in the halo, the efficiency of 


gamma ray production by pp collisions in it is low. (In ad dition, a very long con 


accelerated protons in the halo is needed ~ 10^° years, see 


inement of 


Crocker fc Aharonian 


20111). We 


showed that in this model it is problematic to reproduce the gamma-ray and radio fluxes 
from the FBs, and an additional component of primary electrons is necessary. The magnetic 
held in this model is strongly restricted. The model that reproduces the observed gamma 
and microwave emission from the FBs if the magnetic held is within the range 2.5 to 7 pG. 


Below we analyse an alternative model of stochastic acceleration of electrons in the 
Bubbles. It is known that relativistic electrons are produced by supernova remnants (SNRs) 
which are distributed in the Galactic Disk. These electrons hll an extended region (about 
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several kpc above the plane) of the Ga 


observations (see e.g., 


Berezinskii et al. 


actic halo as fonnc 


1990 


Strong et al. 


from radio and gamma-ray 


20111) ■ However, dne to inverse 


Compton and synchrotron energy losses, only electrons with relatively low energies can 
penetrate into regions high above the Galactic Plane while high energy electrons have shed 
most of their energy before they reach these regions. The rate of energy losses is described 
as dE/dt = —/ii?^ where the parameter n depends on the density of b ackgronnd photon s 


and the strength of the interstellar magnetic held (for details see e.g.. 


Berezinskii et al. 


990l) . For the model of diffusion propagation of CRs presented, e.g., in 


Ackerman n et ah 


(120121) . the length scale of electron mean path length for energies E = 10^^ eV is less than 1 


kpc. Therefore, in order to produce the observed nonthermal emission from the FBs, these 
electrons should be in-situ re-accelerated there. An advantage of this model in comparison 


with the model of stoc 


lastic accelera tion from the background plasma with the temperature 


2 keV fICheng et al. 


2014. paper I ) is that in the case the initial energy of accelerated 


electron is already high, ~ 1 GeV, and the mechanism of re-acceleration needs to increase 
the energy of electrons penetrated into the Galactic halo by three orders of magnitude only. 

Charged particle in the FBs can be accelerated by scattering by MHD-waves 


'Mertsch & Sarkar 

2011 

) or by interaction with supersonic turbulence ( 

Cheng et al. 


2 OI 2 I) . These processes can be described as diffusion in th e momentum space 


the formal equations of 

this process see e 

A- ' 

foDtvgin 

1985; 

Bvkov & Tontvgin 

1993; 

Berezinskii et al. 

199( 

1). We assume 


Bykov fc Fleishman 


or 


1992 


the turbulence can be supplied by active processes in the Galactic Center (GC) when stars 


are captured by the central superma ssive blackho 
can be released by one capture (see ICheng et al. 


e. Energy as hig 


2006 


2007 


1 as hF = 10^^ — 10^^ erg 


20111) . We notice that even 


more energy can be r eleased in the GC, W ~ 10^° erg, if a giant mo l ecular cloud is captured 


by the black hole (see 


Yang et al. 


2012 


20121 ) 


War die fc Ynsef-Zadeh 


(120141) might have found tracers of the last capture of 10^ Mq of gas which occurred ~ 10 


6.5 
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years ago. Similar conclusion was obtained from UV data by 


Fox et al. 


( 2014 ). They found 


indications on a strong outflow from the GC with the velocity > 900 km s“^ that might due 
to past activity of the GC over the last ~ 2.5 — 4 Myr. This time is comparable with the 
age of the FBs. 


Below, we present our analysis of electron re-acceleration in the halo. In Section [2] we 
use the diffusion model to calculate the flux of relativistic electrons emitted by SNRs in the 
Galactic Disk which reach altitudes of about several kpc and stochastic re-acceleration of 
these electrons there up to energies about 10^^ eV. In Section 3 we calculate the fluxes of 
radio and gamma-ray emissions from the region of re-acceleration and compare these results 
with the data derived from the FB. In Section 4 we analyse the effect of convection transfer 
on spectra of accelerated particles and radiation. Section 5 provides a conclusion. 


2. The number and spectrum of re-accelerated electrons in the diffusion model 


In order to estimate the number of re-accelerated electrons and their spectrum in the 
FBs, the kinetic equation should include processes of particle propagation. The reason 
is that the electron sources are in the Galactic Plane while acceleration processes take 
place high above the Galactic Disk. The kinetic equation for the distribution function of 
electrons, f{r,z,p), in this case is 


-V ■ [D(r, 2 ;, p)V/ - u{r, z)f] + 


1 d 2 

p2 QpP 


dp 

dt 


Vu 


P f 


( 


= Q{p^r)5{z) , 


( 2 ) 


where r is the galactocentric radius, z is the altitude above the Galactic plane, p = E/cis 
the momentum of electrons, u is the velocity of the Galactic wind, D and k are the spatial 
and momentum (stochastic acceleration) diffusion coefficients, c{dp/dt) = dE/dt = —pE'^ 
describes the rate of electron energy losses, and Q describes the spatial distribution of 
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cosmic ray (CR) sources in the Galactic plane {z = 0) and their injection spectrum. All 
parameters of this equation are discussed in Appendix 


To dehne regions of stochastic acceleration or re-acceleration in the halo we used 
a hydrodynamic code to simulate the propagation of energy release d in the GC in an 


exponential atmosphere of the halo. We adopted the code PLUTO flMignone et ah 


200711 


and ran hydrodynamic simulations with cylindrical symmetry. We described in the 


introduction that our ic 
SMBH at the GG (e.g., 


ea on the formati on of the FB is a series of star captures by the 


Gheng et ah 


201lh . We simulate the energy released by the capture 


as an explosion at the GG. For illustrative purpose a typical result is shown in Fig. [T] 
Kinetic energy distribution of the gas is plotted in the hgure to emphasize turbulent regions 
as we are going to discuss stochastic acceleration processes in the bubbles. This is the result 
of 100 captures with each energy release 10^^ erg and the interval between two successive 
captures is 10^ years. The distribution shown is the results at 10^ years. 


From the hgure one can see that the shock propagation is mainly in the direction 
perpendicular to the Galactic Plane and the morphology resembles the FB. A layer of 
highly turbulent region is developed close to the envelope of the bubble. The structure 
is similar to that excited by the Rayleigh-Taylor instability (RTI). The development of 
the RTI at the shock front in an exp onential atmosphere has been studied analytical by 


Baumgartner Sz Breitschwerdt 


(120131). 


Similar effect of the RTI has been seen in SNR shocks ("see e.g. 

Hesl 

ter et ah 

1990) 

Numerical calculations near a SNR shock provided by 

Yang & Liu 

2013 

) showed that 


the stochastic acceleration of electrons near a SNR shock by magnetized turbulence may 
dominate over the shock acceleration because most energy of the magnetic helds may be 
generated via the RTF The total energy density of accelerated electrons in this case is 
of the order of energy density of magnetic fluctuations. We cannot exclude that similar 



























mechanism is effective behind the FB shock as iMertsch fc Sarkarl (120111) assumed 


To facilitate discussion and calculation, we show schematically the region of electron 
re-acceleration in the FB where k 7 ^ 0 in Fig. [2] by dark gray. For calculations we used th e 


model parameters taken from the GALPROP numerical program (lAckermann et ah 


2012 


for details see Appendix |^. The momentum diffusion coefficient is taken in the form, 


k{p) = ap^. 


( 3 ) 


The parameter a can be presented as (see e.g., iBerezinskii et al.lIlQOOl) 

C2 

Q! ~ — 


( 4 ) 


where u is the frequency of particle of scattering by magnetic fluctuations of the wave 
number k is, e.g.. 


5Hikf 


( 5 ) 


Here 5H{k) is the strength of magnetic fluctuation with the wave number k, Hq is the large 
scale magnetic field, and ujh = eHQ/m^c. 

In the phenomenological model, which we investigate below, the goal is to estimate the 
value of a from observational data. 


In the simplest case of pure diffusion propagation of CRs when the convection terms 
are neglected (u = 0 ), the number of electrons reaching the re-acceleration region can be 
calculated from 


V-T>(r,z,p)V/ + 4 ^P^ 

pZ fjp 


dp , .df 


= Q{p,r)S{z) 


( 6 ) 


The boundary conditions on the surface of the re-acceleration region are rather 
questionable. The particle may freely escape fro m the boundary as it is assumed for 


the Galactic halo (see e.g. 


Berezinskii et al 


I 990 I) . On the contrary. 


Yang et al. 


fl2ni2h 
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assumed the particle diffusion across the bub ble edge is s 


‘magnetic walls” at the edge as proposed by 


u ongly suppressed, or there are 


Jones et al 


( 20121) . Both effects prevent 


particle propagation through the bubble surface. 


Below we accept for the boundary conditions of continuity for the particle density 
and flux. From Eq. (E]) we calculated numerically the spectrum of SNR electrons which 
reach the altitude z = 5 kpc when the term of re-acceleration is neglected (k(p) = 0). The 
spectrum is shown in Fig. |3]by the solid line. 


The density of high energy electrons needed for the observed gamma-ray flux is shown 
by the shaded region in Fig. [3l The necessary number of high energy electrons can be 
provided by processes of re-acceleration of SNR electrons in the FBs. The spectrum of SNR 
electrons re-accelerated in the FBs is calculated from Eq. ([B]) when the acceleration term 
is included (k(p) ^ 0). However, with the acceleration in the form of Eq. (|3]) the spectrum 
of accelerated particles is too hard (/(p) oc p~^) that is shown schematically in Fig. |3]by 
the dashed dotted line, i.e., our numerical calculations show that too many high energy 
electrons are produced in the re-acceleration region. 


The spectrum of energetic partic’ 
from the acceleration region (see e.g.. 


es can be steepen ed by processes of particle escape 


Cheng et al. 


20121) . Indeed, the momentum spectrum 


of accelerated particles is power-law, /(p) oc p with the spectral index 5 given by 


. 3 

^=2 + 


Q 7" 

' acc 


n. 


(7) 


where the acceleration time Tacc ~ and escape time is t^sc ~ Here Ar^ 

is the thickness of re-acceleration region and Db is the spatial diffusion coefficient in the 
re-acceleration region given by 


Db{p) = 4uV/(6k(p)) • 


( 8 ) 


Here v is the velocity of turbulent motion. Escape processes make the spectrum steeper 
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and thus decrease the number of emitting electrons. 

The spectrum of re-accelerated electrons in the Bubbles calculated for the model 
parameters derived from the observed FB gamma-ray emission (see Section [3]) is shown in 
Fig. |3]by the thin dotted line. The ratio Tacc/Tesc — 8 is shown in next section. 


3. Gamma-Ray and Radio Emission from the FB 


In order to calculate the spectrum of accelerated particles in the Bubble, fb, from Eq. 
([ 6 ]) we calculated the total distribution function of electrons, /(r, z,p), with the acceleration 
term (k 7 ^ 0). Then we calculate the distribution function, /o(r, z,p) when the acceleration 
term is neglected (k = 0). We define the function the excess of electrons, fb, due to 


the acceleration as, fb = f — fo- The procedure is si milar to the su 
gamma-ray flux from the total Galactic emission (see 


Su et ah 


Dtraction of the 


2010 


Yang et ah. 


201411 . 


'B 


Then the FB gamma-ray intensity in the direction of the galactic coordinates {£, b) is 

=-^ J ds J n{e,r)de J p‘^fbir,p) dp. (9) 

where s{£, b) is th e line of sight in the directio n {£, b), and {dfa/de dp)^^.^ is the Klein-Nishina 


cross-section (see 


Blumenthal fc Gould 


19701). We notice that the re-accelerated electrons 


hll not only the region of acceleration but escape into the surrounding medium of the halo. 
Therefore, the total distance of emission is I = Ar -|- \/ Dt where r is the lifetime of emitting 
electrons. 


Our calculations show that in order to reproduce the FB gamma ray spectrum the 
following model parameters are required: the spectral index of accelerated electrons 6 = 4.8, 
the thickness of the re-acceleration region Ar^ = 3 pc, and the characteristic time of 
acceleration a = 2 x 10 “^^ s“^. 
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The calculated FB gamma-ray spectrum and the data from lAckermaun et ahl (120141) 
are shown in the top and middle rows of the left column panels of Fig. 01 As one can see 
there is no complete agreement between the calculation and the data. 

This model has more serious problems to describe the microwave emission from the 
FBs. The calculated spectrum is shown in Fig. 0] (bottom row, left column). As one can 


see, for the parameters derived from the gamma-ray data, the mo del gives ra c 


steeper than v follows from measurements of Planck (see Fig. 


Ade et al. 


io sp ectrum 


(120131)) 


These shortcomings can be alleviated when CR convection propagation is included. 
Below, we present a model of FB that takes convection into account. 


4. Effect of convection on the re-acceleration process 


The effect of convective transfer (Galactic wind) may also 

3e essential 

in the Galaxy as 

it follows from theoretical treatments (see 

Breitschwerdt et al. 

1991, 

2002 

) as well as from 


interpretation of observational data (see 


Bloemen et al. 


19931 ). 


Analysis o 


Crocker et al. 


the radio emission from the c entral Galactic region provided by 


(129111) and 


Carretti et al 


(120131) might indicate a very strong wind 

there. They estimated the w ind velocity up t o 1100 km s“^. Recent three-dimensional 


hydrodynamic simulations of 


Mon et al. 


(120141) showed the existence of strong winds in the 


FBs caused by the past accretion in Sgr A*. They concluded that the wind is collimated by 
the Central Molecular Zone towards the Galactic poles, i.e., perpendicular to the Galactic 
Plane. Thus, we expect that the effect of wind transfer in the GC may not be negligible. 

If the wind velocity is spatially nonuniform in the Galactic halo, then CRs lose 
their energy by the adiaba tic losses. Indeed, the MHD numerical calculations of 


Breitschwerdt et al. 


(jl99lf) showed that the velocity of Galactic wind increases almost 
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linearly with the altitude 2 ;. 


Bloemen et ah 


fll993l) derived the value of the wind gradient 


from the cosmic ray chemical composition. For the wind velocity in the form u{z) = ?)VqZ 


they estimated the gradient value Vq as, Vq ~ 10 s 


-15 ^-1 


Formal solutions for the one-dimensional kinetic equation with the convection term are 
presented in Appendices B and C. 


To demonstrate this effect we solved Eq. 


convection terms. According to 


Crocker et al 


@ w ith the both re-acceleration and 


(120111) we assume that the wind is blowing 


mainly from the central part of the Galactic disk, i.e., from the Bubble region. Therefore 
we took the wind velocity in the following simplified form 

= 3no^6'(ro - r) (10) 

where 9[r) is the Heavyside function and ro = 3 kpc. 


To define the spectrum of accelerated electrons we estimated the number of electrons 
that can reach the altitude about several kpc (shown in Fig. [3] by the thick dashed-dotted 
line). Then we calculated the spectrum of electrons re-accelerated in the Bubbles. The 
effect of adiabatic losses is shifting the spectrum as a whole into the range of smaller 
energies (see the solution in Appendix O and an example in Fig. |8]) and that makes the 
spectrum of re-accelerated particles flatter than in the diffusion model (see the thin dashed 
line in Fig. 12]). 


In the case of leptonic origin of the FB radio and gamma-ray emission, which are 
generated by synchrotron and inverse Compton energy losses of electrons respectively, 
the necessary strength of magnetic held i n the FBs c an be estimated from the simplihed 


equations for these processes presented in 


Ginzburd (119791) . For a power-law spectrum of 


relativistic electrons, N{E) = K^E the magnetic held strength is 


H ~ 


A CWphCTT 


2a(7e 


(4/36) 


(7e-3)/2 


/47rmci/\ 


2/(7e + l) 


( 11 ) 
















~ 13 - 


which is independent from the spectral parameter K^. and the thickness of radiating region, 
if it is the same for the radio and gamma-ray emission. Here a( 7 e) — 0.1, ax is the Thomson 
cross-section, ly. and are intensities of radio and gamma-ray emission from the FBs, e 
and Wph are the energy and the energy density of background photons in the FBs. 


For the FB radio intensity ly = 0.52 kJy sr ^ at the frequency z/ = 23 GHz and 
produced by scattering on optical or IR photons (see 


gamma-ray intensity R = 4 x 10 ® ph cm ^ s ^ GeV ^ sr ^ at = 10 GeV which is 


Cheng et al. 


201 ll) whose energy 


density in the halo is about w ~ 0.2 eV cm we obtain for the electron spectral index 
7 e = 2 that the magnetic held strength is 


H ~ 5 /iG 


( 12 ) 


This value is of the order of one derived below from more accurate numerical calculations. 


The procedure of calculating the spectrum of accelerated particles fb is the same as 
described in Section [2J 

The calculated spectra of gamma-ray emission at different latitudes and radio from 
the Bubbles are shown in Fig. 0] (right column). The best agreement with the data is 
achieved for vq = 10“^^ s“^. The magnetic held strength is H = 3/iG. The parameters of 
the acceleration are the following: the thickness of the re-acceleration region is Ar^ = 60 pc 
and a = 2 X 10“^^ s“^. 


We notice, however, that as it follows from Eq. fITT]) it is problematic to reproduce 
the gamma-ray spectrum in the leptoni c model if we accept the magnetic held strength in 


the Bubbles to be 15 yuG as derived by 


Carretti et al. 


(120131) . The density of relativistic 


electrons estimated from the radio data is too low in order to generate enough gamma-ray 
photons by inverse Compton in the FBs. 


At low Galactic latitudes the contribution of FBs should decrease especially in the high 
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energy range because of energy losses of electrons. On the other hand, at these latitudes the 
contribution from electrons emitted by SNRs in the Disk increases. This effect is shown in 
Fig. IS] where we present the gamma-ray spectrum in the direction of low (10° <| 6 |< 20°) 
and high (40° <| 6 |< 60°) Galactic latitudes. In this hgure we show the IC component 
of gamma-ray emission produced only by the FBs (solid lines) and the total IC emission 
produced by both FBs and SNR electrons (FB-I-SNR, dashed lines). As one can see, in 
the framework of the model the contribution of the FBs to the total gamma-ray flux is 
signihcant at high enough latitudes. At low latitudes the IC emission from SNR electrons is 
dominant and, if we take also into account the gamma-ray component from proton-proton 
collisions, which is very intensive at low latitudes, we conclude that it is almost impossible 
to subtract the FB component from the total gamma-ray flux in these directions. 


The total power, W supplied by external sources (Fermi acceleratio n) which is needed 


to pr oduce high energy electrons in the FBs, can be estimated from (see 


Chernyshov et ah 


20121 ) 


^ - dP’ ( 13 ) 

0 

where p and £ are the particle momentum and the kinetic energy, f(p) is the particle 
distribution function, and k is the diffusion coefficient of the Fermi acceleration. 


We estimate W from the observed FB gamma-ray and microwave fluxes. We take into 
account all processes of electron energy losses as well as their escape from the Galaxy (see 
Appendix [A] for detail). The power estimated numerically is about hF ~ 2 x 10^® erg s“^. 


This estimate is a lower limit for W, because a part of the energy released in the GC 


is also transformed into acce’ 
showed in 


Cheng et ah 


erated protons and the plasma heating in the halo. As we 


(I 2 OIII) processes of star accretion onto the central black hole can 
provide in average about 10^^ — 10^^ erg s“^. Thus tidal accretion supplies enough energy 


for particle acceleration in the GC. 












15 


Finally, we address the question whether the model is able to reproduce sharp edges 
of the Bubbles. As an example we presented in Fig. [B] the longitudinal distribution of 
gamma-ray intensity. As the geometry of the acceleration depicted in Fig. [2] is quite 
schematic, we do not expect complete coincidence between the calculations and the data. 
Nevertheless, the model reproduce qualitatively the effect of sharp edges and neither 


to 


Yane et ah 

2012 ; 

Jones et ah 

2012 ) 


5. Conclusion 

In summary, we studied a leptonic model for the gamma-ray and microwave emissions 
from the FB. The source of electrons is SNR in the Galactic Disk. The electrons propagate 
into the Galactic Halo and are re-accelerated in the high turbulence regions in the FB which 
located several kpc about the Galactic Disk. We have two goals to achieve. The hrst one is 
to get the density needed for the observed gamma-ray flux from the FB, and the second one 
is to obtain the Planck microwave spectrum from the FB. We presented models of different 
level of sophistications. 

To meet the hrst goal, the high energy electron density should be within the shaded 
area in Fig. |3l We introduce the pure diffusion propagation model where we calculate the 
spectrum of electrons (from SNRs in the disk) penetrating into the region of re-acceleration 
and calculate the spectrum of re-accelerated electrons. The re-acceleration process generates 
very hard spectra Therefore electron density we obtained is too high (the thin 

dashed dotted line above the shaded area in Fig. [3]). We then introduce a fast particle 
escape from the re-acceleration region. For small enough escape time the spectrum (thin 
dotted line in Fig. [3]) passes through the shaded area in Fig. |3l 
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Although we got the density of electrons needed for the FB gamma-ray flux, the 
spectrum is too steep to produce the Planck microwave spectrum from the FB and 

we fail the second goal. To remedy, we then include in the diffusion with escape model 
the wind transport and adiabatic losses. The adiabatic losses shift the spectrum to lower 
energy range. Adjusting the parameters of particle escape and wind in the halo we can 
obtain a flatter spectrum which satisfies both goals (gamma-ray flux and Planck microwave 
spectrum). 

The conclusions of the paper can be itemized as follows: 

• Numerical calculations showed that the energy of SNR electrons penetrating into the 
upper halo region is not high enough to generate there the FB gamma-ray emission 
by the inverse Compton. Therefore, further re-acceleration up to energies about 10^^ 
eV is needed there to generate the gamma-ray flux. 

• Re-acceleration (without convection) generates too steep electron spectra. Therefore 
this model is unable to reproduce correctly the gamma-ray and microwave emissions 
from the FBs. 

• There are indications on an intensive outflow of plasma from the GC. The effect of 
the wind leads to adiabatic losses of CRs. We expect that the adiabatic losses make 
the spectrum of electrons in the acceleration region harder than for the case without 
convection. Our calculations show that the gamma-ray and radio emissions of the 
re-accelerated electrons nicely reproduce the Fermi-LAT and Planck datapoints in 
this case. 

• In the re-acceleration model with convection, the gamma-ray flux produced by the 
FBs is more significant at high Galactic latitudes than at low latitudes, see Fig. El 
At low latitudes, the IG emission from the FBs electrons is dominated by that from 
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Galactic disk SNRs electrons. Alternatively, gamma-rays can be produced by pp 
collisions (hadronic model), but this process is not effective in the Galactic halo where 
the gas and GR densities are low. As a whole, at low latitudes the contribution of 
FBs to gamma-ray flux is subordinate to other processes. We should point out that 
gamma-ray production by hadronic model (e.g., pp collisions) could be distinguished 
from leptonic model (e.g., invers e Gompton process of electrons), because pp collision s 


also produce neutrinos (see e.g.. 


2012 


Tavlor et al 


Grocker fc Aharonian 


2011 


Lunardini fc Razzaque 


201411 . 


• An advantage of the re-acceleration model in comparison with that of acceleration 
from background plasma is that in the hrst case the energy of SNR electrons should 
be increased in the FBs by three orders of magnitude only, while in the second case 
electrons are accelerated from the thermal plasma with the temperature about 2 keV, 
i.e., nine order of magnitude increase is needed. 


• We compared the efficiency of electron acceleration for the two models: acceleration 
from background plasma presented in Paper I and the re-acceleration presented in 
this work. To do this we included both processes into the kinetic equation (E]) for 
electrons. Our numerical calculations showed that for the same momentum diffusion 
(acceleration term) and the time of escape the re-acceleration mechanism is more 
effective for the production of high energy electrons in the FBs (see Fig. ED- 
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A. Parameters of the CR kinetic equation 

The distribution function of electrons, f{r,z,p), is derived from the kinetic equation 
(see also Eq. (|2])) 


-V- [D{r,z,p)Vf -u{r,z)f] + 


1 cl 2 

p2 QpP 


dp V-u \ df 

— - p t — K(r, Z, p)-:— 

dt 3 V ^dp 


= Q{p,r)S{z) 


(Al) 


Relation between momentum and energy losses is 

dp 1 dE 
dt 13 dt ' 


(A2) 


where (3 = v/cis dimensionle ss velocity of the electron and the rate of electron energy losses 


can be presented as (see e.g.. 


Ginzburg 


197911 : 


— - f—'] f—'] 

dt \dt ) \dt ) \dt \dt ) 


(A3) 


where the rate of Coulomb losses in a neutral medium and in a fully ionized plasma are 
respectively 

( ^ ^ ^ 10“^^ GeV s~^ X nHl3~^ [log(7 “ 1)(7^ ~ 1) + 20.5] , (A4) 

\ dt J cn 
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) = ’^•62 X 10 GeV s ^ x uhii^ ^ [log( 7 /n^f//) + 73.6] 


(AS) 


Here 7 = E/{mc^) = 1/\/l — /9^ is the Lorenz-factor, uh and uhii are densities of neutral 
and ionized hydrogen in the disk and in the halo. 

The bremsstrahlung losses in the interstellar medium are 

I —— ) = 5.1 X 10“^® GeV s“^ X (uh + tihiA'J , (A 6 ) 

\dt Jbr 

and the synchrotron and inverse-Gompton losses (in the Gompton limit) can be presented 


as 


=1.6x10 GeV s ^ X -I-7^/3^ 


(A7) 


The Klein- Nishina cross-section for inve rse Gompton scattering of high energy electrons can 


be found in 


Blumenthal fc Gouldl ( 1970 ). 


Ackerman n et al. 


We take the values from the GALPROP code (see 
parameters in Eq. flAip . 

The distribution of magnetic held in the Galaxy was taken in the from 


H = 6 fiG X exp 


2 OI 2 I) for other 


5 kpc 20 kpc 


(AS) 


as the densitv of t 

le interstellar 

f Ackermann et al. 

(2012 

)• 


The distribution of SNRs in the disk was taken in the form 

1.2 


Q{p,r) = Q{p) X 


r 


exp ( -3.22- 


r 


8 kpc / \ 8 kpc 

As follows from the analyses of radio and gamma-ray th e spectra of CR e 


protons have a depletion a t energies below 1 GeV (see e.g., 


2012 


Dermer et al. 


Strong et al. 


2011 


(A9) 
ectrons and 


Neronov et al. 


20131 ). Here we take the source spectrum of electrons as a broken 
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power-law from 


Ackermann et a1. 


(1201211 


-3.6 


Q{p) oc < 


P 


-4 4 

p , 


p 


-6 


if pc < 2.2 GeV 

if 2.2 GeV < pc < 2.1 TeV 

if pc > 2.1 TeV 


(AlO) 


The spatial diffusion coefficients in the Galaxy is assumed to be a function of the 
momentum only, 

/ DC \ *^-33 

Dip) (^) . (All) 

where Dq = 9 x 10^® cm^ s—1 for 8-kpc halo. 

Parameters inside the acceleration region were derived from the FB gamma-ray and 
radio emission. The momentum diffusion coefficient has the form 


/ \ 2 


(A12) 


where acceleration rate Ob is estimated from calculations, and the spatial diffusion coefficient 
there is 

4y2p2 


Db{p) = 


6Kb{p) ’ 

where v is the ve l ocity of turbulent motions which provide stochastic acceleration. In 


(A13) 


Ackerman n et al. 


(120121) this acceleration is due to particle interaction with Alfvenic waves. 


then V = va = H/ \/47rp. 


The boundary conditions were taken in the form 


9//9r = 0, at r = 0 

df /dz = 0 , at z = 0 (A14) 

/ = 0 , at the boundary of the Galactic halo. 

At the boundary of the acceleration region, continuity of the particle function and flux are 
enforced. 
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B. Re-acceleration in a divergent wind flow 


In this appendix we present a solution to a one dimensional re-acceleration problem in 
divergent flow. The governing equation is 


df ^ df 9 f ^df\ f du\ pdf 1 d f 2 df 


dt ^dz dz\^dz) \dz} ?> dp p"^ dp \ ^dp} ^ 


(Bl) 


We consider a speciflc case where u = Uqz/H, D = Dq, k = kqp'^ = a‘^p‘^/9DQ, Q = —f/r, 
where Uq, Dq, kq, a and r are constants. H is the characteristic length of the system or 
flow. We seek steady state solution with boundary conditions: f = s{p) at z = 0 and / = 0 
as 2 ; —)■ cxo. Introducing dimensionless quantities 


Mo = 


Huo 

Dn 


~ Ha . Dqt z ^ p [3 dp 3 

^ = T7> P = —) h=/-—= -logp, B2 

Do H po J a p a 


and 


71 

v) = v) exp — {a^ + Mo)J = v) exp(-z/r/), 


the steady state of Eq. flBip becomes 


ft, f’ 


(B3) 


(B4) 


and the boundary conditions for g become: g = s{p)p^’'H = S{p)e'''^ at ^ = 0 and 5 ^ = 0 as 
f —)■ cxo. We can solve Eq. flB4p by Fourier transform with respect to p. We obtain an ODE 

+ + = (B5) 

subject to boundary conditions: g = S{u) + iu) at ^ = 0 and ^ = 0 as ^ 00 . Here the 

Fourier transform pairs are 

/•OO poo 


and 


9ii,v)= du , g{^,uj) = — g{^,p) e dp . 

J — 00 J —00 


/ OO 1 POO 

S(u) e‘“'' doj , S(n) = — I S(r,) e-*™ <ir,, 


(B 6 ) 


(B7) 
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The solution of g{^,u) can be written in terms of parabolic cylinder function 




2-’^/2 n-n 


S{u + iu)exp ( ) Vni^/^0 > 


where r(a;) is the gamma function, is the parabolic cylinder function, and 

1 f Un Mn 1 

n = -—{^ + — + ^ + - + u 
Uq \ 2 4 4(7^ r 


Dehne 




2-n/2 _ /I 


r 


exp 


UoC" 




and the Fourier transform pair 


/ OO 1 poo 


The solution can then be written as 

1 


9ii,v) = 


Siv')e‘'^'g{^,r,-r,')dr,', 


2tt 


and 


/(e, v) = 7^ S{v') V - i) dll'. 


27T 


(B8) 

(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 


C. Analytical Solution of the One-Dimensional Wind Equation 


The one-dimensional equation for relativistic electrons can be presented in the form 
fp‘ ]\i f) f) 

- DoE'^-^ + 3^0^ {zN) - — {fiE^ + voE) N = QoE-'^^Siz). (Cl) 

Here DqE’^ is the coefficient of electron diffusion in z direction, V{z) = ?)VqZ is the wind 
velocity in z direction, dE/dt = —fiE"^ is the rate of synchrotron and inverse Compton 
energy losses, QoE~z°S{z) is the source function of electrons in the Galactic disk. 


Introducing variables 

VqE Zd 



rj = 


z 


T = 


Zd = 


r“(l -|- ry ’ 


(l + r)3’ 


(C2) 
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and the function 


K = 


ZdVo / no 


70 


(1 + ^r 


N. 


Qo 

we obtain the standard one dimensional diffusion equation 




(C3) 


(C4) 


where F{t) = r 


— T-7o+a-2 


(1 + tY with r expressed as a function of t b y inverting 


t = f dr/ [r“(l + tY]. The Green function of Eq. flC4p can be found in 

imi- 


Morse fc Feshbach 


The solution of this equation can be obtained in the analytical form 


Qoh E-(t'o+“/2) rdrorY-" 

N = -TT / -- exp 


VDoVq (1 + fiE/voY 


Vo 


EY D,E»{l + fiE/voY 


(C5) 


where 


S = 


dx 


x°‘{fiE/vo + xY 


(C6) 


TO 


The effect of adiabatic losses is a shift of the spectrum into the region of low energies. As 
an example we presented in Fig. [8] the spectrum of electrons at the altitude 2 ; = 6.6 kpc 
calculated for the velocity gradients Vq = 2 x 10“^® s“^ (upper line) and Vq = 5 x 10~^^ s~^ 


(botto m line). The production spectrum of electrons by SNRs was taken from 


(120 llh as 


Strong et ah 


Q{E) = ( 0(4 GeV - B) + g(E - 4 GeV). 


(C7) 
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Multiple Explosions 100xl0''(53) erg 
8 5 



r(l<pc) 


Fig. 1.— Distribution of kinetic energy of gas from the hydrodynamic simulations of Fermi 
Bubbles. In this simulation, we model the energy release by a star capture by SMBH at GC 
as an explosion at GC. The hgure shows the case of 100 captures with each energy release 
1053 interval between two successive captures is 10^ years. The distribution is 

the result at 10^ years after the first capture. Kinetic energy is plotted to emphasize the 
turbulent regions. A layer of strong turbulence occurs close to the Bubble boundary. 
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r 


Fig. 2.— The schematic picture of the Galactic halo (light gray) and the re-acceleration 
region of the FB (dark gray). In the model we put h = 8 kpc, = 3 kpc, Zb = 5 kpc and 
Zq = 3 kpc. The thickness of re-acceleration region, Ar^ is estimated from the results of 
numerical calculations. 
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Fig. 3.— The spectrum of SNR electrons accelerated in the Bubbles at the altitude 5 kpc for 
the diffusion coefficient D = 10^® cm^ s“^, /r = 2 x 10“^® s“^ GeV“^ shown by the thick solid 
line when convection term is neglected. The same spectrum of SNR electrons penetrating 
into the re-acceleratio n region for the convect ion model with the velocity gradient vo = 10“^® 


s ^ (as derived e.g. in 


Bloemen et ah 


fjl993[B is shown by the thick dashed-dotted line. The 


re-accelerated spectrum of electrons in the Bubbles without escape term (oc E~^) is shown 
schematically by the thin dashed-dotted line. The same for the case, when escape processes 
from the re-acceleration region of the thickness Ar^, = 3 pc are taken into account, is shown 
by the thin dotted line. The spectrum of re-accelerated electron when convection transfer is 
essential is shown by the thin dashed line. The density of electron needed for the observed 
gamma-ray flux from the Bubbles is shown by the gray region. 
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E (GeV) E (GeV) 






Fig. 4.— Spectra of gamma-ray and microwave emission produced by re-accelerated elec¬ 
trons in the FBs calculated for the diffusion model (left column) and for the case when the 
effect of convection is included (right column). The top and middle rows are gamma-ray 
emission for the latitude range 20° < 40° and 40° < 60°, respectively. The bottom row 


presents the spec 


Ackermann et al. 


rum of microwave emission 


(120141) for gamma and from 


rom the Bubbles . Datapoints were taken from 


Ade et al 


(120131) for the microwave spectrum. 
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Fig. 5.— Spectra of IC gamma-ray emission in the direction of low (left panel) and high 
(right panel) Galactic latitudes. The spectrum of total IC emission (by FBs and SNRs of 
the Galactic Disk) is shown by the dashed line, and the contribution by FBs only is shown 
by the solid line. 



Fig. 6.— 
40° < & < 


Longitudinal distribution of the FB gamma-ray in t ensity for the latitude range 


50°. Datapoints were taken from 


Ackerman u et ah 


fl2014h . 
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Fig. 7.— Comparison of acceleration efficiency for the cases of acceleration from a back¬ 
ground plasma and for the case of re-acceleration for the same parameters of momentum 
diffusion and electron escape. 







(snunqjB)N 
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E(GeV> 


Fig. 8.— This figure is an illustration of adiabatic loss in divergent flows. Spectra of electron 
at the altitude 2 ; = 6.6 kpc for the velocity gradients vq = 2 x 10“^^ s“^ (upper line) and 
Uo = 5 X s“^ (bottom line). 
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